D2 layers on MgO(0 0 1): Simulation study

In response to recent helium atom scattering (HAS) and neutron scattering results, Monte Carlo simulations and perturbation theory calculations have been performed for D₂ on MgO(0 0 1). Monte Carlo simulations predict that D₂ molecules form a series of interesting structures, p(2 × 2) → p(4 × 2) → p(6 × 2), with coverages Θ = 0.5, 0.75, 0.83 respectively, and followed by a formation of a top layer of p(6 × 2) unit cell symmetry. The three types of mono-layers are stable up to 13 K, whereas the top layer still exists up to 10 K. This is in partial agreement with the neutron scattering and HAS results that report c(2 × 2), c(4 × 2) and c(6 × 2); they agree in terms of coverage and stability, but disagree in terms of symmetry. A quantum mechanical examination of the D₂ molecules’ rotational motion shows the molecular axes are azimuthally delocalized and hence the simulated structures are c-type rather than p-type. These calculations also indicate that ortho-D₂ and helicoptering para-D₂ prefer cationic sites, while cartwheeling para-D₂ prefers anionic sites.     

Phase behavior and fluid-solid surface tension of argon in slit pores and carbon nanotubes

Phase behaviors of argon in several types of cylindrical and slit pores are examined by grand-canonical Monte Carlo simulations. Condensation processes in single- and multi-walled carbon nanotubes along with those in hard-wall tubes are compared. Effects of the pore size on pressure-tensor components, the fluid-wall surface tension, and the adsorption are also compared for the different fluid-pore interactions. The chemical potential at which the fluid begins to condense in the single-walled nanotube is greater than that in the multi-walled nanotube by an amount nearly equal to the difference in the potential-well depth of the fluid-pore interaction, and the adsorption isotherms overlap each other almost completely for narrow pores and partially for wider pores. Similar analyses are performed for slit pores of two different hydrocarbon models.     

System-parameter dependence of the metallic phase of the non-doped 2D Hubbard model

Naito et al. reported that some non-doped T′-214-type compounds drive high-T_c superconductivity. The compounds are considered to be metallic since on-site Coulomb energy U is moderate and the Fermi surface is much deformed in these compounds. In order to confirm this picture and extract electronic structure information, we have examined the phase diagram of the metallic state of the 2D Hubbard model as a function of U and t′ (with t″ we fixed at − t′/2 here; t′ and t″ are the second- and third-neighbor transfer energies, respectively) by means of the variational Monte–Carlo method. We employed a Jastrow-type Gutzwiller trial wave function. In the studied range of U = 2–12, the boundary value for |t′| at which SDW disappears increases almost linearly with U. Jump-wise transition to the Mott insulator state was not observed. Using the boundary curve and experimental band parameter values, we estimate U  5 for T′-214 compounds. Preceding works are discussed in the last part.     

Momentum dependence of pseudogap and superconducting gap in variation theory

To consider the origin of a pseudogap and a superconducting (SC) gap found in the high-T_c cuprates, we evaluated the momentum dependence of the singlet gap corresponding to the pseudogap and the SC gap in the t–J model, using an optimization variational Monte Carlo (VMC) method. In the underdoped regime, the singlet gap is significantly modified from the simple d_x²-y²(d)-wave gap (∝ cos k_x − cos k_y) by the contribution of long-range pairings. Its angular dependence at the quasi Fermi surface is qualitatively consistent with those experimentally observed in both hole and electron-doped cuprates. On the other hand, a SC gap is almost unchanged, preserving the original simple d-wave form. Thus, it seems that the incoherent part of the singlet gap mainly influences the forms of observed gaps.     

Study on CYCIAE-100 radiation field and residual radioactivity

The accelerators should be properly designed to make the radiation field produced by beam loss satisfy the dose limits. The radiation field for high intensity H- cyclotron includes prompt radiation and residual radiation field. The induced radioactivity in accelerator components is the dominant source of occupational radiation exposure if the accelerator is well shielded. The source of radiation is the beam loss when cyclotron is operating. In this paper, the radiation field for CYCIAE-100 is calculated using Monte Carlo method and the radioactive contamination near stripping foil is studied. A method to reduce the dose equivalent rate of maintenance staff is also given.     

Monte Carlo simulation of the property of a scintillation bar in the multi-neutron correlation spectrometer

To perform a kinematically complete measurement of the dissociation reaction for neutron-rich nuclei, a multi-neutron correlation spectrometer is proposed at Peking University.A Monte Carlo simulation code based on GEANT4 is developed for a single scintillation bar which processes not only the energy deposition but also the light propagation in the scintillator and the light collection and conversion to signal at the end of the bar in a realistic way. The simulating method is described in detail in this paper, and the timing and position resolutions and detector efficiency are studied based on the simulation and compared with the experimental results.A new method of crosstalk rejection has been demonstrated to be important for the design of the whole spectrometer.     

Total reaction and fusion cross sections at sub- and near-barrier energies for the system 7Li + 28Si

Fusion cross sections are extracted for the ⁷Li$ + $²⁸Si system, via reaction cross section and transfer measurements at sub- and near-barrier energies ( E _lab = 5.7 to 14MeV). The energy evolution of transfer to reaction cross section ratios is determined with the aid of CDCC calculations, which subsequently allows the deduction of fusion cross sections at sub- and near-barrier energies. It is shown that fusion can be well represented in a BPM context. Fusion cross sections are compared for the systems ⁷Li$ + $²⁸Si and ⁶Li$ + $²⁸Si, the latter studied previously, and are found to exhibit different strengths. Last, the direct channels determined at 13MeV, are found to be dominated by a 2n -transfer mechanism.     

Fluctuation-theory constraint for extensive entropy in Monte-Carlo simulations

The entropy per particle in most Monte-Carlo simulations is size dependent due to correlated energy fluctuations. Guided by nanothermodynamics, we find a constraint for the Ising model that enhances the fluctuations and lowers the free energy, while making the entropy homogeneous, additive, and extensive. Although the average interaction energy becomes size dependent, the resulting distribution of energies provides a mechanism for the heterogeneity found in the dynamics of many materials.     

Transitions between secondary structures in isolated polyalanines

Monte Carlo simulations of gas-phase polyalanine peptides have been carried out with the Amber ff96 force field. A low-temperature structural transition takes place between the α-helix stable conformation and β-sheet structures, followed by the unfolding phase change. The transition state ensembles connecting the helix and sheet conformations are investigated by sampling the energy landscape along specific geometric order parameters as putative reaction coordinates, namely the electric dipole μ, the end-to-end distance d, and the gyration radius R_g. By performing series of shooting trajectories, the committor probabilities and their distributions are obtained, revealing that only the electric dipole provides a satisfactory transition coordinate for the α↔β interconversion. The nucleus at the transition is found to have a high helical content.     

晶粒棱长、尺寸与拓扑学特征之间关系的Monte Carlo仿真研究

采用Potts模型Monte Carlo方法研究了晶粒棱长、尺寸与拓扑学特征之间的统计关系.结果表明,晶粒棱长与晶粒面数之间呈线性统计关系,并且平均N面体晶粒模型和Poisson-Voronoi组织均支持该结论.不同时刻的晶粒长大仿真数据表明,在准稳态晶粒长大阶段晶粒棱长的分布具有自相似性.个体晶粒的平均棱长随晶粒面数（或晶粒尺寸）的增加而逐渐增大,这说明一些理论模型中采用的“不同面数的晶粒平均棱长均相等”的假设具有局限性.仿真数据和纯铁实验数据均表明,晶粒尺寸与晶粒面数之间的统计关系表现为一条单调递增的凸曲线. 关键词： 晶粒棱长 晶粒尺寸 拓扑学 Monte Carlo仿真